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^ ' After giving an overview of recently invented methods for simulating lattice QCD 

, at small li/T, we discuss some results for bulk thermodynamic quantities of QCD 

i matter coming from those methods. We focus on the transition line and the critical 

f"*^ ■ endpoint in the QCD phase diagram. 
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1. Introduction 



Having a profound and quantitative knowledge of the QCD phase diagram 
and equation of state is of great interest for meany different topics: A few 
pi 1 micro seconds after the big bang, the early universe went through the QCD 

O . phase transition; Every heigh energy experiment has a certain trajectory 

in the phase diagram, which may also cross the transition line; Finally the 
cores of cold neutron and quark stars may be in a color superconducting 
phase, also part of the QCD phase diagram. 



$_i ' During the last two decades lattice simulations have proven to be 

the only quantitative approach to the QCD phase transition. At vanish- 
ing baryon chemical potential ((1b — 0), lattice simulations provided de- 
tailed informations of several bulk thermodynamic quantities, for instance 
the transition temperature (T c ps 173 MeV) and critical energy density 
(e c fa 0.7 GeV/fm 3 ) 1 . The sign problem, however, has disabled simulations 
at non-zero chemical potential. For /i > the determinant of the fcrmion 
matrix (detM) becomes complex. Standard Monte Carlo techniques using 
importance sampling are thus no longer applicable when calculating observ- 
ables in the grand canonical ensemble according to the partition function 

Z(ji) = f X>f7 detMQx) exp{-S G }. (1) 

Note that all physical observables are real and therefore also the imaginary 
part of detM = |detM| exp{i6*} vanishes in average. What influences the 
physics and courses the problems are the fluctuations of the complex phase 
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9. Several new techniques have been developed to overcome this problem 
in the region of small /z/T, where the fluctuations are rather mild. For a 
recent overview see for instance Ref. 2. 



2. Methods for small fi/T 

2.1. The Taylor expansion method 

One conceptionally very simple idea is to calculate observables (O) at it > 
from a Taylor expansion around /i = 0: 

O(fi) = c + c 1 fi+^c 2 fi 2 + --- (2) 

Since on the lattice all quantities are given in units of the lattice spacing 
(a), the expansion parameter is jj, = a/j, = N 1 r l (^i/T). Here N t is the 
number of lattice points in the temporal direction. The idea goes back 
to the calculation of the quark number susceptibilities in Ref. 3. By this 
method the response of meson masses 4 as well as the pressure and further 
bulk thermodynamic quantities 5 ' 6 have been studied. Besides derivatives of 
the observable itself, the calculation of derivatives of In detM with respect 
to /i are required, the first two nontrivial coefficients in Eq. (2) are given 

by 

/dO\ / dlndetMV 

/ d 2 0\ /<90<91ndetM\ / <9 2 lndetM\ / \ / <9 2 In detM \ 

The derivatives have to be taken at /to = 0. Note that due to a symmetry of 
the partition function (Z(fi) = Z(—u)) all odd coefficients in Eq. (2) vanish 
identically. For the same reason we have (dlndetM/9/t) = at /to = 0. 
We explicitly use this property in Eq. (3) to derive the coefficients. The 
advantages of this method are: 

(i) Expectations values only have to be evaluated at fi = (no sign 
problem) . 

(ii) All derivatives of the fermion determinant can be expressed in terms 
of traces by using the identity In detM = Tr In M. This enables the 
stochastic calculation of the coefficients by the random noise method, 
which is much faster than the evaluation of the determinant. 

(iii) continuum and infinite volume extrapolations arc well defined on a 
coefficient by coefficient basis. 
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On the other hand it is not a priori clear for how large /J./T the method 
works and how large the truncation errors are. Furthermore one is strictly 
limited by phase transitions, since phase transitions are connected with 
discontinuities or divergencies. An estimation of the convergence radius of 
the series gives a lower bound on the applicability range and thus also a 
lower bound to the phase transition line in the (T, //) plane. 



2.2. The reweighting method 

In principle one can also calculate the exact expectation value (O) (fi) 
from an ensemble generated at fi = 0. This is possible by giving ev- 
ery configuration an additional weight according to the reweighting factor 
R(fi) = dctM(/t)/dctM(0). We have the identity 

<O)(/OH0(£W)>/W)>. (4) 

where the expectation values on the right hand side are to be evaluated at 
ft = 0. This method goes back to Ref. 7. For large reweighting distances, 
however, the generated ensemble at ft, = shares less and less configurations 
with the target ensemble. This is the so called overlap problem. The overlap 
can be improved by the multi-parameter reweighting technique 8 , where the 
reweighting factor is a function of fi and the lattice coupling fi. In the 
(/3, fi) plane one can reweight along lines where the overlap measure 9 (a) is 
maximal. One of those lines is the transition line. In Ref. 9 the generated 
ensemble at (/3o,0) was defined to have the overlap a with the ensemble 
at ((3,fi), if a fraction a/2 of the ensemble gives the weight 1 — a/2 to 
the expectation value at (/3,/t). In Fig. 1(a) the volume dependence of 
the overlap measure along the transition line is shown. In Ref. 9 it was 
empirically found, that the half width fi X j 2 of a scales with fi X j 2 oc V~~' ', 
with 7 « 1/3. 

To reduce the numerical effort, which is connected with the calculation 
of the reweighting factor (R), one can expand R in a Taylor series around 
/to = as it was done in Ref. 10. In lowest order we have 

DM 



In R = TrM 

ofi 



xA + C(A 2 ) (5) 



In this expansion all odd terms are purely imaginary, whereas all even terms 
are real 10 . The complex phase of the determinant is thus given by 



9 = 6 + 0(fi 3 ), with 9 = ImTrM- 1 ^ 

O/J, 



x fi. (6) 

Ao=0 
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Figure 1. (a) The overlap measure along the transition line for different volumes 9 . 
Shown lattice sizes are 6 3 X 4, 8 3 X 4, 10 3 X 4 and 12 3 X 4 (from top to bottom), (b) The 
approximated phase S vs. the exact phase 8 (Rcf. 11). 



Again the problem of the determinant calculations is reduced to calculations 
of traces. In Ref. 11 it was empirically shown that and coincide quite 
well for /i<0.3 (see Fig. 1(b)). In the region /t>0.3, where one also faces 
the overlap problem, the lowest order expansion seems to break down. 



2.3. Analytical continuations 

At imaginary chemical potentials, the fermion determinant is real and pos- 
itive, thus simulations by standard Monte Carlo techniques are possible. 
Results on the imaginary jlj axis can be analytically continued to the real 
fi R axis. It is especially easy to convert a Taylor series in fij, expanded 
around jlo = 0, into a Taylor series in jln- Since the series has only even 
powers of /t, due to the the symmetry Z(fi) = Z(—fi), one only has to switch 
the sign of every second coefficient (c 2 — > —c 2 ,cq — > — c 6 , . . . ). There is how- 
ever another symmetry of the partition function which limits the analytic 
continuation. Due to the periodicity 12 Z(ptR, Hi) = Z(ur, /i/ + 27rT/3) sim- 
ulations with /ij > will only give excess to the physical region /j,r<itT/3. 
This corresponds to ^s<500 McV. This method was also used to map 
out the phase transition line 13-15 . One should note, that for this method 
neither an evaluation of the determinant, nor any of its derivatives are re- 
quired. In order to determine a Taylor series of an observable one needs 
however many different simulation points for several values of fii, to in- 
terpolate by a certain Ansatz. In future, one should think of combining 
methods 2.1 and 2.3 to get a better extrapolation along the fiR axis. 
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2.4. The density of states method 

An alternative to the importance sampling technique used in most Monte 
Carlo simulations is the density of states method. Here one reorders the 
path integral of the partition function in the following way: first expectation 
values with a constrained parameter will be calculated. I.e. one exposed 
parameter ((f)) is fixed. Expectation values according to the usual grand 
canonical partition function (Zqc) can then be recovered by the integral 



where the density of states (p) is given by the constrained partition function: 



With ( we denote the expectation value with respect to the constrained 
partition function. In addition, the product of the weight functions /, g 
has to equal the correct measure of Zqc- fg = detMexp{ — Sg}- This idea 
of reordering the partition functions is rather old and was used for gauge 
theories 16 and QED with dynamical fermions 17 . For QCD the parameter 
(f> is normally set to be the plaquette 18 : <j> = P. In Ref. 19, however, 
the DOS method was constructed for the complex phase (</> = 9). Within 
the random matrix model, the authors of Ref. 20 used the quark number 
density ((f) = n q ). 

The advantages of this additional integration becomes clear, when 
choosing </> = P and g(U) = 1. In this case p(<j)) is independent of all 
simulation parameters. The observable can be calculated as a function of 
all values of the lattice coupling (3. If one has stored all eigenvalues of the 
fermion matrix for all configurations, the observable can also be calculated 
as a function of quark mass (m) and number of flavors 18 (Nf). 

Note that this method does not solve the sign problem. It is, however 
supposed to solve the overlap problem. Nevertheless, it is also possible to 
combine the DOS method, with the reweighting method 2.2, by reweighting 
the constrained expectation values in the case of g(U) ^ 1. For large 
reweighting distances an overlap problem is introduced then once again. 

3. The critical temperature 

On the lattice, the critical temperature (T c (p)) can be calculated from 
the critical coupling (f3 c (jl)). We define the function (3 c (jl) as the line in 
the (/3,/t) plane, where susceptibilities peak. Together with the lattice 
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Figure 2. (a) The transition line T c (fi) for diferent number of flavors and masses in 
the range m/T = 0.1 — 0.4, for the light quark mass, (b) The mass dependence of the 
transition line, for Nj = 3 (Ref. 23). 



beta function ((3(a)) one can calculate Taylor coefficients of the critical 
temperature, for the first nontrivial coefficient we have 

dT c i dp c (fi) f opy 1 

d( M 2 ) iV 2 T c (0) 9(A 2 ) Vda] ■ W 

This quantity was calculated by several groups 10, 13 - 15 > 21 . 22 j with methods 
2.1-2.3. Although different quark masses in the range of m/T = 0.1 - 0.4, 
lattice sizes and actions has been used, the results are in good agreement 
(see Fig. 2(a)). In first order, the transition line for Nf = 2 and m/T = 0.4 
is given by 10 

T c (^ q )/T C (Q) = 1 - 0.070(35) (fi q /T c (0)) 2 . (10) 

As one can anticipate from Fig. 2(a), the mass dependence of T c (fi) is 
rather small. Comparing Nf = 3 results 23 (Fig. 2(b)) for m/T = 0.4 and 
m/T — 0.02 one finds however a significant mass dependents. As expected 
the curvature of the transition line becomes stronger with increasing mass. 
Using the perturbative asymptotic /3-function in Eq. 9 gives 

TM f 1 - 0-025(6) (^)',m/T = 0.4 _ ^ 
T <=(°) |l -0.114(46) (j^y) 2 , m/T= 0.02 

We note that an additional quark mass dependence is hidden here in the 
transition temperature at \i q — 0, i.e. T c (0), which is used to normalize 
the transition temperature at \i q ^ 0. Furthermore, taking into account 
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violations of asymptotic scaling in the /3-function will lead to a further 
increase of the curvature of T c (p q ). 

This increase of the curvature can be understood as a result of much 
larger phase fluctuations for the smaller quark mass. A detailed analysis 
of the phase fluctuations can be found in Ref. 24. To make this more 
clear we have also calculated the transition line as a function of a vector 
chemical potential (jj,v = p u = — Md)- This corresponds to reweighting 
with the absolute value of the determinant, i.e. the influence of the phase 
is dropped. In this case the mass dependence is much less significant: 

Tc fa v ) _ Jl -0.026(5) (fifc) ,m/T = 0.4 



r c (o) 



,2 (12) 
1 - 0.033(6) (t^j) ,m/T= 0.02 



Since for the larger quark mass the two transition lines are in good agree- 
ment, we have T c (p v ) > T c (p q ) for m/T = 0.02. This result seems to be 
a bit counter intuitive, since we know that the onset transition of isospin 
matter at T = takes place at fiy = rather than at \i q = m^/3. 

However, this result can also be confirmed by the DOS method. In Fig. 3 
we presents results 25 from a 4 4 lattice and Nf = 4. Simulations have been 
performed at finite isospin density with \iy = 0.3, according to Ref. 26. 
The DOS (p) was defined according to Eq. 8 with g = |detM| exp{— Sq} 
and <j> = P. In Fig. 3(a) we show p, the phase suppression factor {cos(9)) 
and their product as a function of the constrained Plaquctte P. We find, 
that the phase fluctuations in the low temperature phase (P<2.8) are much 
larger than in the high temperature phase (P>2.8). From this we can ar- 
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gue, that including the phase factor, the critical temperature will be lower 
than without the phase factor. Indeed, after calculating the susceptibilities 
from those distributions as a function of /3, we find a shift in (3 C towards 
lower temperatures (see Fig. 3(b)). Moreover the value of (i c is in good 
agreement with an earlier result form the reweighting method 8 a . 

4. The critical end-point 

The determination of the critical end-point on the transition line, separat- 
ing a first order region form a crossover region is a difficult task. At present 
there are two different strategies used for this calculation. An infinite vol- 
ume extrapolation of Lee- Yang Zeros 21 ' 22 and a finite size scaling analysis 
of susceptibilities and Binder Cumulants 27 ' 28 ' 13 . At this stage the most re- 
liable calculation for the critical endpoint and physical quark masses comes 
form Rcf. 22. The result is \i B = 360(40) MeV and T = 162(2) MeV. Note 
that those values are, however not yet free from corrections coming from 
the finite volume and lattice spacing. 

In the case of 3 degenerate quarks the critical mass (m c (0)), where the 
critical point sits at /i = 0, was determined. In this quantity huge cut-off 
effects were found. It was determined for standard staggered fermions 27 
and pA improved fermions 23 , both on N t = 4 lattices. The outcomes for 
the critical pseudo scalar mass are m£ rit = 190(20) MeV, 67(18) MeV, re- 
spectively. Binder Cumulant method also admits the determination of the 
universality class of this critical point, which was found to be that of a 
3d-Ising model 27 . 

In addition the mass dependence of the critical point turned out to be 
rather large. A correct calculation of the critical point may thus only be 
possible, if the quark masses are known very precisely. The first Taylor 
coefficient of m c {p) was calculated in Ref. 13. The result of a combined fit 
for several lattice sizes and imaginary chemical potentials is 



where m c (n) is the critical quark mass. This coefficient does not look 
unnaturally small, since it was expanded in the natural units (fi/nT), it 
will, however give in turn a rather large first coefficient for the critical 
chemical potential as a function of the quark mass {jjb c {m q )). 



a Agreement was found after a linear extrapolation A — * 0, with A being a small parameter 
in front of an isospin symmetry braking term as defined in Rcf. 26 




(13) 
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Figure 4. (a) The absolute values of the Taylor coefficients of the pressure 6 , (b) The 
convergence radius of the pressure 6 . 



Another interesting method to extract information about the critical 
end-point is the convergence radius of any bulk thermodynamic observable 
at p = 0. As mentioned in Sec. 2.1, the Taylor expansion is limited by a 
phase transition. Thus for any T < T cn d po int = Tb the convergence radius 
(p) should remain finite, whereas for T > Te it should be infinite. Since 
physical observables are even in p, we have 

p = lim p n = lim y/\c n /c n+2 \. (14) 

n— >oc n— >oc 

In Fig. 4(a) we show the absolute value of the Taylor coefficients of the 
pressure 6 . The results are for Nf = 2 and m/T = 0.4. In Fig. 4(b) the first 
two terms in the series (Eq. 14) are drawn in a (T, /^-diagram, in order 
to compare with the phase transition line. As one can see in Fig. 4, for 
values of T/T c (0)>0.8 the radius seems to increase. This corresponds to a 
critical chemical potential of p E w (3 - 3.6)T w (600 - 700)MeV. This is 
a reasonably large value for m q /T = 0.4. 

Note that the pressure and quark number density, which can be con- 
structed from the coefficients in Fig. 4(a) is in good agreement with result 
from method 2.2 (Rcf. 9). Moreover the pressure is like the transition line, 
almost perfectly described by the leading order in p/T. For a further dis- 
cussion of the pressure from Taylor expansion see Ref. 29. 



5. Conclusions 

For small p/T several methods have been proposed and tested. The agree- 
ment between different methods is good. The phase diagram for small p/T 
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is slowly evolving and one can even address the question of locating the 
critical point. All methods are bases, however, on ^ = simulations and 
thus extrapolate. Moreover they give no solution for the sign problem and 
the results are not yet extrapolated to the continuum and infinite volume. 
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